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In this review article we give a brief overview on some recent progress in quark pairings in dense 
quark/nuclear matter mostly developed in the past five years. We focus on following aspects in 
particular: the BCS-BEC crossover in the CSC phase, the baryon formation and dissociation in 
dense quark/nuclear matter, the Ginzburg-Landau theory for three-flavor dense matter with Ua{1) 
anomaly, and the collective and Nambu-Goldstone modes for the spin-one CSC. 



I. INTRODUCTION 

T-H ^ It is well known that fermion pairing is the underlying mechanism for superconductivity and superfluidity. In a 
■ fermionic system, the weak and attractive interaction between two fermions leads to the formation of Cooper pairs 
' at low temperatures, which is well described by the Bardeen-Cooper-Schrieffer (BCS) theory. In strong interaction 
<^ |. quarks and gluons are elementary particles which are described by quantum chromodynamics (QCD). There are also 
^ ' quark pairings called color superconductivity (CSC) first proposed in cold dense quark matter by Barrois Bailin 
and Love [2], and further developed by many others (for reviews, see }3l-[l3j). However color superconducting quark 
matter may not be relevant in collisional experiments in the foreseeable future. The reason is that this exotic phase 
of matter requires extremely high baryonic densities and relatively low temperatures. In nature, such conditions 
may be realized in the cores of cold c omp act stars, i.e. in relatively old remnants of supernova explosions (for an 
introduction of compact star, see e.g. |14|). Among some major developments o n q uark pairings in the past decade 
are comprehensive studies of important CSC phases on the phase diagram |15l 420[ : the color-flavor-locking (CFL) 
state J2l|. 2^. the two- flavor (2SC) state [23 - [28| . the single flavor pairing state |27l. [29l - [43t . mis-matched pairings 



12, liJlland many other phases (Zsl. . The construction of a theory for the CSC at weak couplings is another 
important progress [21413, 113, [MiJ^ 

In this review article we will give a brief overview on some recent progress in quark pairings in dense quark/nuclear 
matter mostly developed in the past five years. We will focus on following aspects in particular: the BCS-BEC 
" ^' crossover in the CSC phase, the baryon formation and dissociation in dense quark/nuclear matter, the Ginzburg- 
Landau theory for three-flavor dense matter with J7a(1) anomaly, and the collective and Nambu-Goldstone modes for 
— . ^ the spin-one CSC. Random phase approximation and Dyson-Schwinger equation are used to obtain the propagating 
, modes of diquark pairings. The BCS-BEC crossover are investigated within the boson-fermion model and the NJL- 
type model. The baryon formation and dissociation in different phases are studied by regarding a baryon as composed 
, of a quark and a diquark in the NJL-type model. With a nonlocal extended NJL model one can obtain the constituent 
' quark mass which is momentum dependent. The bubble diagrams are automatically convergent, providing an effective 
^-H ' confinement mechanism. In a three-flavor NJL model including the axial anomaly, a low temperature critical point 
^-H is found due to the coupling between chiral and diquark condensates. The collective modes in the spin-1 CSC are 
^ ^ analyzed in the Ginzburg-Landau approach. 
. 

H II. BEC-BCS CROSSOVER IN A RELATIVISTIC BOSON-FERMION MODEL 

The Cooper pairs are formed by two fermions at the Fermi surface in weakly attractive channel, therefore they 
have spatial extension called the coherence length which is much larger than the mean inter-particle distance. In 
a strongly coupling regime, the Cooper pairs are bound to bosonic molecules and condense in the ground state to 
form the so-called Bose-Einstein condensation (BEG). The Fermi surface disappears and no fermion degree of freedom 
remains in this situation. Although the features of BEG and BCS are very different, there is no phase transition but 
just a crossover between them. Recently many experimental advances have been made in cold atom system. With the 
help of an applied magnetic field, the effective attractive interaction between the atoms can be tuned via Feshbach 
resonances and the BCS-BEC crossover can be observed. 

The BEC-BCS crossover in the CSC was first studied by Nishida and Abuki [55*, '56'| with the NJL model and 
followed by others [57l - l65| . Inspired by those previously used in the context of cold fermionic atoms, one can directly 
extend the boson-fermion model to a relativistic version. In this case both the fermion and difermion channels are 
included as fundamental degrees of freedom with total fermion number density fixed. Tuning the effective interaction 
strength the crossover between the BEG and BCS can be studied [57, 64]. In a non- relativistic system, people normally 
use the scattering length to characterize the strength of the attractive interaction between fermions. In the relativistic 
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boson-fermion model, a crossover parameter x is defined by the difi'erence between the square of efi^ective boson mass 
and boson chemical potential. The inverse of x play the similar role of scattering length. A varying magnetic field 
can also lead to the relativistic BEC-BCS crossover [6a|. Starting from any initial state at zero field, with a ultrahigh 
magnetic field the system always settles into a pure BCS regime. 

For the NJL type model for the BCS-BEC crossover in the CSC phase (55l . [56l . [ssl . [soj . since the fundamental 
degrees of freedom are fermions, the Cooper pairs are introduced with the help of bosonization procedures like 
Hubbard-Stratonovich transformation. Some nonperturbative tools such as random phase approximation and Dyson- 
Schwinger equation are needed to determine the propagating modes of the boson field. By taking the phase shift to 
the non-relativistic limit (p — > 0), an effective scattering length between fermions can be derived in relativistic case. 

We review the boson-fermion model for the BCS-BEC crossover with the total fermion number density fixed (57l.[6^. 
The Lagrangian respects the global U{1) symmetry, which consists of free fermion and boson parts, Cf and and 
a Yukawa interaction part Ci, 



C = Cf+Cb + Ci, (1) 



with 



= giipipcil^tp + ^*ipi7^'>Pc)- (2) 

The fermions and bosons degrees of freedom are described by the spinor tp and the complex scalar boson field (p. The 
charge conjugate spinors are defined by ipc = Cip^ and V^c = ^r^C with C = 17^7°. The fermion/boson mass is 
denoted by m/mf,. The Lagrangian is invariant under the U{1)b transformation ?/; — e"*"')/', Lp — >■ e^'"(/9. Considering 
a system with the total U{1)b charge conservation, the boson chemical potential is chosen to be twice the fermion 
one, Hh = 2/z. Therefore, the system is in chemical equilibrium with respect to the conversion of two fermions into one 
boson and vice verse. This allows one to model the transition from weakly-coupled Cooper pairs made of two fermions 
into a molecular di-fermionic bound state, described as a boson. The bosonic field are separated into two parts as 
if — if^ + (j)^ where ip^ is the expectation value of the bosonic field in vacuum or the zero mode and cf) is the nonzero 
mode. The condensation is denoted as A = 2g(^o- The fermionic field can be rewritten in the Nambu-Gorkov (NG) 
basis, ^ = {tp^ipc)'^ , ^ = {'^','ipc)- The bosonic field can also be rewritten as {4>R,(j)i)^ , where and </)/ are the 
real and imaginary part of the complex bosonic field respectively. Then the Lagrangian can be cast into the following 
form 

C = \^S-^^ + 4>i)D-' {\])+ {^n, 4>i)^ ( + ■ (3) 



Here we have re-arranged the boson-fermion interaction in the third term. The inverse tree level fermionic and 

bosonic propagators are denoted as S*^^ and -D^^, {Tr., f is the boson-fermion vertex with Tr = iy/2gj^ai'^ and 

Tj — —i^^g^^a^'^ , where the a^^ ^-re the pauli matrices in the NG basis. In the following the gap A is treated as 
real for simplicity. 

With the Lagrangian density one can work out the effective potential in two levels: mean field approximation 
and the two particle irreducible (2PI) approach |f64l|. Li the mean field approximation, only the condensate 
contribution of the bosonic field is included while the fluctuation of the field is neglected. In the 2PI approach the 
calculation is done in the Cornwall- Jackiw-Tomboulis (CJT) formalism [67|, in which all 2PI diagrams are counted so 
that the fluctuation of the bosonic fleld can be calculated self-consistently. The effective potential in the mean field 
approximation and in the 2PI approach read 

Tmf = -/(A) + iTrlni^-i - ^Trln^-^, (4) 

fap/ = -/(A) + i{Trln2?-i+Tr(2?-ip-l) 

-Trln5-i - "Ix^S-^S - 1) - IV^Piip.S)], (5) 
he condensate 

which in the present case is 



where /(A) — ''^"s — is the condensate contribution. V2P1 includes all 2PI contributions to the effective potential. 



« -iTr{P,r.iTr[f,5f,5]}, (6) 
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where Tij are the boson-fermion vertex, i,j = R,I correspond to the real and imaginary components of the bosonic 
field, T) and S are the dressed bosonic and fermionic propagators derived by the Dyson-Schwinger equations (DSE) 
for bosons and fermions. In the C JT formalism the DSE are given by taking derivatives of the effective potential with 
respect to propagators, 

^"^'^ 0, %i = 0, (7) 



which lead to the DSE, 



SV ' SS 



, ST2P1 



V-^ = D-^ -2^^ = D-^ +U{V.S), (8) 
oV 

^ 5-i+2^^^^5'-i-I](I?,5), (9) 
06 

where the self-energies for bosons and fermions in the CJT formalism are given by 

n^Trlf.^f,^], E(P,5)=Tr[P,f.5f,]. (10) 

The Rainbow-Ladder truncation is applied to the calculation for self- energies, i.e. replacing the dressed boson-fermion 
vertices and propagators by the bare ones. Substituting the DSE into the 2PI effective potential, 

f2P7 = -/(A) + i{Trln[Z?-i(l + I?n)]-Trln[5-i(l-5E)]-Tr(I]5)} 

« -I{A) + ^{TrlnD-'^ -TrlnS-^ +TiiDn{V,S)]}. (11) 

In the last step an expansion in terms of self-energy has been made leading to a partial loss of self-consistence but 
reducing the numerical complexity. From the last line of Eq. (Ilip . by comparing the mean- field with 2PI effective 
potential, one can find the beyond- mean- field contribution to the effective potential is included in -—T2Pi. 
From the DSE for the fermion, Eq. ([9]), the 11-component in the NG basis reads 

- [{S-')ll - Sll - (5-l)i2((5-l)22 - ^22)-\S-%i] 

« [{S-')n~^ii-iS-')i2i{S-')22)-\S-%i]. (12) 

In the last line E22 is neglected since it is sandwiched by (5^^)i2 and (5^^)21 which are proportional to A. Sn is 
the 11-component of quark self-energy En — {TT:[DijTiSTj]}ii. Then we have 



-Sii -{S-')i2{{S-%2)-\S-') 



21 



Po + go - - eEp+q ^ ^ g Po - A* - eSp 
+ A^ 

E ^i /Ar- (13) 



The assumption q ^ p is used in the second step. = \/q'^ + is the boson energy in vacuum. Apg is the 
pseudo-gap defined as 

d^q 1 



= 16g' 



(2^)3(po+/ib)2-(i^,T 

d^q l + lB{El-iJ.^) + fB{E\ + iit,) 



(27r)3 2Eb^ 

2 / d\ fB{E^,~f^,) + fBiE',+f^,) 
^ ' (27r)3 2Eb ■ ^ ' 
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(a) 



D(«>). 



(b) 



Figure 1: Diagrams for the gap and pseudo-gap. The indices 1 and 2 are for NG space. 



In the last step, the pseudo-gap are renormahzed by directly removing the divergence part 1/{2E^) which conies 
from vacuum. From the fermionic self-energy Eq. (|13p . the lowest order contribution of the fluctuation is to add 
a pseudo-gap term to A'^. Here we define the pseudo-gap as the correction to the fermion self-energy at the 
static limit, see Fig. in In the region near and above the superconductivity critical temperature Tc, the di- fermion 
fluctuation will contribute to the fermionic spectral density near Fermi surface by bringing in two bumps structure, 
but it is not a real gap since the fermionic excitation is not forbidden between the bumps. It has analytical structure 
from which one can compute the density of states from the di-fermion fluctuation. In the region below Tc both the 
gap and pseudo-gap contribute. Here the pseudo-gap effects are approximated by A^^ for numerical simplicity. The 
dressed fermion propagator is also approximated by replacing A^ with A^ -I- A^ 



From Eq. (jlOp . the bosonic self-energy have four components corresponding to z,j — I,R. There are simple 
relations among them: Hrr/ii = Hq ± Hi, Hir = Uri = 0, where Hq — —8g'^Tr['y^ Suj'^ S22] is the self-energy in the 
normal phase (A = 0), and Hi — — 8(7^Tr[7^52i7^52i] is proportional to A^. The momentum integrals are implied in 
Tr. Also we have Drr = Djj for the bare bosonic propagator. Inserting the self-energy to Eq. ^ we obtain r2p/ as 

Tap/ = -^TT[DRRURR]~^Tr[Djjnjj] = -TT{DRRUo). (15) 
Then the effective potential in the mean field approximation and the 2PI approach are given by 



fMP = -E/ ^[4 + 2THl + e-^^/-)] 



4e 



45 



2^^ J i2n) 



^2PI = Tmf— r2P/, (17) 



2rin[l - e"'^^/^], (16) 



where = Vfc^ + ^ti^ — e/i and — y^(^^J^~f~A^ are fermionic excitation energies in normal and condensed phases 
respectively. In the following we use to denote the thermal dynamic potential to replace the effective potential f . 
The total charge density n = is fixed to a constant. The fermion number density for fermions, condensed/excited 
bosons, and the 2PI component are 



2fiAl 

d^k 



PbO 

Pb 



m = — , (18 

n oil 



which satisfies 



Pf + PbO + Pb + Pr2 = 1- 



(19) 
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Figure 2: Left panel: the chemical potential fi (blue dotted line) and the gap A (red dashed line) at zero temperature as 
functions of x. Right panel: the condensed boson and fermion fractions with x. The units for fj,Q, Aq and T are GeV. 



The gap equation is settled with the saddle point condition of the free energy density F = + /in, that is 

dF _ dn dn dti dn _ dn _ 

dA~ dA^'dildA^'^dA' dA^^' ^^^^ 

The gap A, the chemical potential fi and the pseudo-gap Apg can be solved simultaneously from the density and the 
gap equations (|19l20p and the pseudo-gap equation (jl4p . In the mean field approximation T2P1 and Apg are set to 
zero and the number of the equations is reduced to two. 

In the mean field approximation, the renormalized boson mass mtr can be obtained as 

-.---4, ^\A.T.,.o^ml-4g^J (^^p^- (21) 

In the CJT formalism to calculate the renormalized bosonic mass nibr one should work with the bosonic DSE ^ for 
self-consistency. Then the pole equation in the homogeneous limit with gap set to zero is given by 

det2?-i|p=A=A„=o - dct[i?-i +n(po,p)]|p=A=A,,=o = 0, (22) 
with which the dressed boson mass mi,r is determined as 

rribr = {po + fJ-b)^ ^ ml + Reno{po - i^), (23) 



-(Po+A^b)ry-Imno(po-«^) =0, (24) 

where po is the pole position and 77 is the width of the dressed boson propagator determined by the imaginary part 
of the pole equation ([^^ . The condition 77 = defines the bosonic dissociation boundary (T*,/i*). By solving the 
real and imaginary parts of the pole equation simultaneously, the dissociation boundary turns out to be very simple 
nibr = 2™. In the present model the {T*,ii*) are functions of the bare boson mass mb- From Eq. ([3]), with a fixed 
boson- fermion coupling constant g, mt, is equivalent to fix the coupling constant in a pure fermionic model. Assuming 
the boson is stable, then we get r/ = and Eq. becomes 

mlr ^ ml +no{po), (25) 

_ 2 

nibr together with the bosonic chemical potential serve as the crossover parameter x = _ I!i^_^jiik , xhe parameter x 
can be varied from negative values with large modulus (BCS) to large positive values (EEC). In between, xq is the 
unitary limit. Therefore, x behaves as the scattering length. 

In Fig. [51 the calculation is done in mean field approximation for A, /i and bosonic/fermionic charge fractions as 
functions of x. The parameters are set to T = GeV, m = 0.2 GeV, g — 4 and A = 1 GeV. xq = J i^^J'y^ ^k'^+m'' 
an upper limit of x which ensures non-negative bosonic occupation number. In the right panel of Fig. [21 the charge 
fraction of the thermal boson is always equal to zero since the temperature is zero. From the negative value to positive 
value of X, the system goes through a crossover from the BCS side to the BEG side. 
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Figure 3: The fermionic chemical potential /i (blue dashed line), the gap A (red solid line) and pseudo-gap (black dotted line) 
as functions of mi,. Right panel: fractions of fermions, condensed/thermal bosons and 2PI component. The units for nib A*, A 
and A„„ are GeV. 



At high temperatures, both the condensed and thermal bosons should be considered. The results for the calculation 
up to 2PI are shown in Fig. [3l The parameters are chosen to be T = 0.14 GeV, m — 0.28 GeV, g — 1.8 and A = 1 
GeV. With a small boson bare mass m;, the system is in a strong coupling regime where the fermion chemical potential 
/i is low and the condensed bosons are dominant. When m;, goes larger, both the gap A and pseudo-gap Apg decrease 
while the chemical potential /i increases. The pseudo-gap Apg is small due to a low temperature chosen here. The 
fraction of condensed bosons becomes smaller and the fermionic degree of freedom becomes more important. From 
the left panel of Fig. [31 the 2PI contribution is dominant in the large mt, region indicating the strong interaction 
between fermions and bosons. 

The fluctuations may change the superconducting phase transition to be of first order, such as the intrinsic fluctu- 
ating magnetic field in normal superconductors ^68i] . or the gauge field fluctuations in color superconductors [54], due 
to the fact that fluctuations bring a cubic term of the condensate to the effective potential making the Landau theory 
of continuous phase transition invalid. In the present case, a nonzero cubic term of A is generated in the bosonic 
fluctuation r2p/, leading to a first order phase transition. 



It is well known that an applied magnetic field can tune the BEC-BCS crossover in cold atom system, because the 
magnetic field can adjust the effective interaction between fermions via Feshbach resonance. In a relativistic fermion 
system, by applying a magnetic field with a magnitude near the energy scale of the interaction, the properties of 
the system will also be affected. The magnetic fields on the surface of pulsars are about 10^^ ^ 10^'^ G , and for 
magnetars they are about 10^'* ~ 10^^ G |69l - [7]| . In the core the magnetic field can be even stronger. In heavy ion 
collision experiments at RHIC and LHC, the background magnetic field generated in non-central collisions can be 
about 10^* - 10^^ G at RHIC/LHC energies In the future low-energy experiments at RHIC, NICA and FAIR 

which is targeted to probe dense and cold nuclear matter, the magnetic field are also expected to be very strong. Such 
a strong magnetic may bring some significant effects to the matter, for example, the chiral magnetic effect |72l. [t^. 
The strong magnetic fields have significant effects on the quark pairings in the CSC phase |75l - [78| . Recently the 
magnetic field tuning of the BEC-BCS crossover has been studied |66j. 

One can extend the boson-fermion model discussed above to study a system with oppositely charged fermions 
= (V'ii'02) and neutral scalar bosons [66]. An external magnetic field is also included in the model and coupled 
with the fermions. Then the fermion part and interaction part of Eq. ([1]) are 



where q denotes the charge of the fermion. is the vector potential of the external magnetic field. As discussed in 
Sec. ini the Lagrangian is invariant under the U(1)b transformation — >■ = e~*"\l/, ^ tp' = e*^"</3. Hence the 
simple relation /it, — 2/i ensured by chemical equilibrium remains. In order to describe the BEG of these molecules, 
we also separate the zero- mode of the boson field ip and replace it by its expectation value (p = (tp) , which represents 



III. RELATIVISTIC BCS-BEC CROSSOVER IN MAGNETIC FIELD 



Ci = ¥'*c(«7'^cr2)* + ¥'**(«7^cr2)^'c, 



(26) 
(27) 
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the electrically neutral difermion condensate. The mean-field efi^ective action is then 



- I |2 -ml I ^ |2, (28) 

where the fermion inverse propagators of the Nambu-Gorkov positive and negative charged fields 5'+ — (?/'2, V'lc)^ 
and ^ ^ = ipic)'^ a-re given by 



with 

[Gf^^o]''{x,y) = [zT^^i) - m ± mtI^'I^^ - y), (30) 

and n^^^ = id^ ± qA^. Without loss of generality, the magnetic field can be chosen along the z-axis with A'^ = 
(0, Bxi, 0, 0). By using Ritus' transformation to momentum space, the effective potential at zero temperature reads 



" = E E^w / ^^^3^^+ 4 +4;^E / ^^P^^^ 

e=±lfc=0 •'° e=±l-'0 



(31) 



in which d{k) = (1 — ^) denote the spin degeneracy of the Landau levels, fc = 0, 1, 2... is the Landau level (LL). The 
energy dispersion of fermions and bosons is given by 



ee(fc) = ymT^^ij^pT^, e = ±l (32) 
We = ^Jp^ + ml- 2e/^, e = ±l, (33) 



respectively. = -^p^ + 2\q\Bk + is the energy of free fermions in the magnetic field. The parameters of the 
model are the momentum cutoff and the fermion mass. The total fermion number density can always be chosen to 
be at X = where the fermion number fractions of fermions and bosons are equal. Here we choose: the momentum 
cutoff is chosen to be a Gaussian type exp[— (p| -|-2|q|i?fc)/A^] with A = IGeV, and the fermion mass m = 0.2A. With 
the effective potential, the gap equation 951/9 A and the density equation —dn/dfi = n can be simultaneously 
solved. 

In this model there are two variable parameters: the bare boson mass rub and the magnetic field B. From Fig. [21 
with B = and mb is small the system is in BEG regime, a large rrib corresponds to BCS regime. Varying mf,, the 
system makes crossovers between these two regimes. In the following the boson mass irib is fixed since the effect of 
magnetic field can be shown. 

Fig. [1] shows the chemical potential, the gap and fermion number fraction of the condensate bosons and fermions 
as functions of B with nib — 0.8 GeV [6gj. As the magnetic field is weak on the left end of the panels, the system is in 
BEG regime with the condensate fraction much larger than the fermion fraction. Increasing B till h\{qB / m?) ^ 0.1, a 
de Haas van Alphen behavior will be present. But in this case there is only one oscillation with small amplitude on the 
curves for /i and A, since the system is in the BEG regime in weak magnetic field with the gap about 73 MeV, while 
the de Haas van Alphen oscillation favors small B and the amplitude is suppressed by a large gap. By setting a large 
enough value for the bare boson mass to let the system start in BCS side, by tuning the magnetic field the oscillation 
will be more obvious and one can find several crossovers in the fermion number fraction. Gontinuously increasing the 
magnetic field, a pure BCS state settles down on the right end of the panels. This process is a BCS-BEC crossover by 
varying the magnetic field, but the origin is different from that in non-relativistic case where magnetic field is used 
to tune the Feshbach resonance. 

The mechanism of the crossover can be explained by the energy dispersion of the fermion shown in Fig. [S] [6^. 
In the left panel of FiglSJ the LLs with fc < 3 contribute to the BCS component on which the minima is located at 
P3 7^ 0, while the one with k — 'i contributes to the BEG one. Since the fermion energy splitting between different 
LLs and the density of states of each LL are all proportional to V ei3, when the field increases not only the energy 
levels become more separated, as seen from the figure, they can also accommodate a larger number of particles. As a 
consequence, when the field increases, the number of occupied LLs reduces, or in other words, the magnetic field will 
press the fermions to lower LLs. Hence, with increasing the field the contribution from the lowest LL becomes more 
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Figure 4: The fermionic chemical potential and the gap (upper panel), and fermion number fractions (lower panel) as functions 
of ln(g_B/m'^). Starting from the BEC side where the BEC component is much larger than the BCS component at small 
magnetic field on the left. The system crossover to a pure BCS state at large magnetic fields on the right. 
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Figure 5: Positive energy component of fermion dispersion relation with fixed mt and B. The LLs fc = 0, 1, 2, 3. 



important and finally dominant. The higher levels shown in the middle and lower panel of Fig. [S]are likely not to 
contribute in strong fields. That is the reason why in a strong magnetic field the system is in the BCS regime. 

If the relativistic BCS-BEC theories discussed in the literature have any relevance for the physics of neutron stars 
and the future low-energy heavy ion collision experiments, we should study the effects of the magnetic field on the 
BCS-BEC crossover within a more realistic model, as extremely strong magnetic fields are expected to be present. 



IV. DIQUARK PROPERTIES AND BCS-BEC CROSSOVER IN DENSE QUARK MATTER 



The interaction between two quarks in the anti-triplet channel in color space is attractive, leading to Cooper pairs 
of quarks at extremely high densities and low temperatures. Due to asymptotic freedom of QCD, the interaction is 
weak and the Cooper pair wave function has a correlation length that exceeds the inter-particle distance. However, 
as the density is lowered, the interaction strength increases and the Cooper pair becomes more localized. Eventually, 
Cooper pairs will form tightly bound molecular diquark states, then the diquark BEC regime is settled. Since the 
interaction between quarks is strong, it is argued that the diquark fluctuation is large around the critical temperature 
Tj, and a pseudo-gap may be observed. When the temperature rises, the pseudo-gap will becomes smaller and finally 
disappear as the dissociation of diquark fluctuation. Recently there are some work that focus on those issues within 
the NJL model. 

One considers a 2SC case in the NJL model. The Lagrangian reads 

Cnjl = V5(^7M9^-"^o + Ai)V' + G,[(Vi^/')^ + (^/5^75^V')^] 

+ Gd [i'il^T2Jatpc] [V'ci7^T-2 JaV-], (34) 

where tp &iid ipc a-re quark field and its charge conjugate respectively. T2 and t are the Pauli matrices in flavor 
space, and {Ja)bc = ~i^abc{oL = 1,2,3) denote the antisymmetric color matrices. Gs/Gu are coupling constants for 
quark- anti-quark/quark-quark channel, which together with the momentum cutoff A and the bare quark mass TOq are 
the input parameters of the model. The hat of mo means the bare quark mass in flavor space. 
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To introduce the chiral condensate and the diquark degrees of freedom, one can use the mean field approxima- 
tion which is equivalent to the Hubbard-Stratonovich transformation in non-relativistic case. The diquark field 
can be decomposed into the condensate part Aa = 2Gd {i>il^T2Ja'>Pc) and the fluctuation part then one gets 
'ipij^T2Ja'ipc = ^2°G • phase the superconducting gap are chosen as Ai 2 — and A3 7^ without loss of 

generality. Then the Lagrangian ((34|) becomes 



Cnjl « -i'I'5-i*--^^|Aj2_G',(a„+ad)^ 



2 4Gd 

-^i^lR + vlI) + l^^a^Ta^■9. (35) 

The quark fields can be expressed in the NG basis, = (i/jjipc)'^ and ^ — {ij;,ipc)- <^u,d are chiral condensates and 
can also be derived in the mean field approach. The inverse propagator then reads 

{ i7'^2 JaA P^Y' -^i-f^~m J' ^"""^ 

where the quark mass in flavor space is to = (too + to^) • 1 / with rriq the chiral condensate —2Gs{<Tu + a-d)- In present 
case only the scalar quark-quark channel is considered. The original complex diquark field has been decomposed into 
two real bosonic fields: the real part cpaR and the imaginary part (fai with ipa = -^(VaR + 'i^ai)- The last term in Eq. 
P5p is a Yukawa type quark-quark-diquark vertex with the index i — I, R, with which the diquark dynamic properties 
can be studied. TaB.i — "^7^''2>/a'''i^2^i where t(^2^ are the Pauli matrices in the NG space. There are additional two 
tadpole terms: (/9qA* and (p*A, but in a self-consistent theory all the tadpole terms should cancel themselves. One 
can prove that they are canceled by the tadpole terms of the one-loop diagram generated by the terms '4)i^^T2Jai'cf*a 
and (paipci^^T2Ja^P, 



loop 



^A* + l/^Tr[.,,r2J.5.] 



4Gu 4 Jk 

= 0. (37) 

This relation is satisfied due to the gap equation in the mean field. Similarly the tadpole terms corresponding to cp* 
can be proved to cancel each other. For clarity the tadpole terms are not included in the Lagrangian. The full diquark 
propagators are derived via the Dyson-Schwinger type equation, 

ATabo,p) = -^-n,a(po,p), (38) 

where po = i^mrT are the Mastubara frequencies (n — 0, ±1, ±2...) and 11.; ^ are the diquark self-energies, which have 
the following properties in the 2SC phase: n^y/.a = ^(ng ± HJ), and = 7^ Hg, = SasUl The expression 
for the Hq I are 



n;,'' = -g' f TT[S22iK)-f\2Jl.2SlliP + K)j^T2Jl,2] 

(27r)3^ 2eier Po - e[et ~ ^1 



7, 1 TP e |Cfc,p+fe, 



= -g2 [ TT[S22{Kh''T2J3Sn{P + K)-f\2.h] 

Jk 
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Figure 6: The phase diagram obtained within the two flavor NJL model. 



f Tr[5l2(i^)7S Jl.25l2(P + Kh^T2Jl,2] 
J K 



0, 



(41) 



5' / Tr[5l2(i^)7V2J35l2(P + i^)7'T2J3] 
J K 



' K 



^3 



(27r)3 eie'ie^e-^ 



(42) 



where the quasi-quark energies are = eEk — /i, e| = y^^|)^ + A^ with = ^1 + m^. Ck,p+k is the product of 



,k-(p+k)H 



'-], where summation is imphed over e, e', ei, e'j^. 



energy projectors, Ck,p+k = [1 - 

With the full propagator of diquarks, one can study the stability and dissociation properties of diquarks. In Fig. [51 
the red solid line separates the chiral symmetry broken phase from the symmetric phase (indicated by xSB/xSR)] 
CSC denotes the color-superconducting phase. In the chiral symmetry broken phase and CSC phase, with a strong 
diquark coupling Gd, due to nonzero mass gap rrig and the CSC gap A^, stable diquark poles can be found in the 
window (— 2mg,2mg) and (— 2Aci,2Aci) respectively in p = limit. In x^B phase (with a large Gd) and normal 
phase, there is a boundary below which diquark pole equations + Ili,a{po,p) — have solutions. That line 
is defined as the diquark dissociation boundary, see the blue dashed lines for three values of the diquark coupling 
constant, Gd = 3.11, 3.8, 4.025 (in units of GeV~^) in Fig. [6l The corresponding regions in Fig. [6] are filled with light 
blue, green, and magenta color, respectively. These poles also exist in the CSC phases, however, for the sake of clarity 
there is no color in the CSC regions. The CSC phases for Gd = 3.11, 3.8, 4.025 are bounded by the red solid line and 
the dash-dotted lines from bottom to top, respectively. Note that the diquark coupling constants we have chosen here 
are in the weak-coupling or BCS regime. As we increase Gd, Bose-Einstein condensation of diquarks can take place 
in the region below the dissociation lines, provided the bare quark mass is nonzero ^551. IstI IssI. Isil. [6^ [6^. FtqI. Isoj . Note 
that in Ref. [6l| , a vanishing decay width was imposed as an additional criterion for the location of the dissociation 
boundary. 

In Fig. [71 the diquark spectral densities are presented for four points: A {T ^ 0.03, /i — 0.25 GeV), B {T = 0.03, 
fj. = 0.33 GeV), C {T ^ 0.03, = 0.36 GeV) and D [T ^ 0.15,/i = 0.36) in the phase diagram Fig. [B] The spectral 
density for diquarks is defined as following 



Pi,a(w,p) 



1 



(43) 



In each panel there are three curves corresponding to three momenta |p| — 0,0.2,0.4 GeV. The diquark coupling 
constant is chosen as a weak one Gd = 3.11GeV-2. The upper three panels of Fig. [71 correspond to points A, B, 
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Figure 7: Spectral densities for diflFerent values of T and p,. The upper panels correspond to (T, /iq) =(0.03,0.25), (0.03,0.33), 
(0.15,0.36) GeV, respectively. The lower panels are for (T, /Xq) = (0.03,0.36) GeV (inside CSC phase). The indices in the 
spectral density refer to I,R and colors. 



D, from left to right respectively. All the six components of spectral density with indices i = R,I and a = 1,2,3 
are identical outside CSC phase. At point A (in X'^B phase), with a weak diquark coupling, there are only broad 
bumps above two times quark mass gap (|a;| > 2mq). The point B is located below the diquark dissociation boundary 
in normal phase, where diquarks have poles corresponding to stable diquark resonances. The point D is above the 
dissociation boundary in the normal phase, and no diquark pole exists. The broad bumps in the diquark spectral 
density indicate unstable diquark resonances. The lower three panels are for ifj and ipa diquark field with color 
indices 1,2,3 at point C in CSC region. The /, R diquark field with color index 3 are gapped while the others are not. 
From the expressions of diquark self-energy, five Nambu-Goldstone modes are recognized, those are the I/R fields 
with color indices 1 and 2, and the / field with the color index 3. They can be directly proved with the fact that the 
diquark pole is located at (a;,p) =0 GeV due to the gap equations for A^. In the lower-middle (spectral density for 
the (fii^s field ) and lower-right panel (spectral density for the </5//fl_i/2 field), one can find the five Nambu-Goldstone 
modes at (a;,p) = GeV. The lower-left panel is for ipR,3 field without Nambu-Goldstone mode. 

To evaluate the thermal diquark contribution to the thermodynamic potential, Abuki has developed a method to 
express the potential in terms of diquark spectral density [s^ . The total thermodynamic potential is decomposed into 
two terms — Qmf + ^fiuc, where flMF is the mean field potential and ^fiuc is the thermal diquark contribution. 
The full diquark propagator at some coupling constant 4Gd — Q is obtained with the dispersion relation 



pL(w,p) 



(44) 



The thermodynamic potential is given by 



flue 



■4Gr 



dp 


log - 


(27r)3 


dg f 


dp 


J 


(27r)3 


dg f 


dp 


G'J 


(27r)3 



duj 



(r = /i = part) 



OJ — lUJri 



(T = /i = part) 



(45) 



The advantage of the propagator in terms of the dispersion relation is that the summation over Matsubara frequency 
can be analytically performed. The fluctuation part can be decomposed into two pieces: the Nozieres-Schmitt-Rink 
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term 



2 J (27r) 
and the quantum fluctuation term 



^NSR = ^§2 j / c^w^p,^^(w,p)[ni3(w) +6i(-w)] 

/oo 
du;^(5i,„(w,p)[nB(a;) + ^(-a;)], (46) 
-OO A „ 



1 f dp 



1 dg f dp z""" , V- A G . 
= -2/0 g^J J^J_J-J:^pU-'P 

I f dp f°° duj^ , 
= -2/ (2^LV5^'^"'^"'^ 



sgn(a;) 



(47) 



in which the integral over the spectral density is represented by the phase shift deflned as 

'45^ -ni,a(w + i77,p) 



/ 

^0 



dg g i. 
Pi,a = o 



2' 



4^ -ni,„(a;-zr?,p) 



5i,a(w,p). (48) 



Here Apf^ = pf^ — (T = /z = part) and AiJj^a = (5j,a — (T = /z = part), ns is the bosonic distribution function. 
^NSR is the contribution from the thermal fluctuation which vanishes at zero temperature, while Slqfi remains at zero 
T. In further calculation the quantum fluctuation contribution is neglected, and only the thermal effect is considered. 
The charge conjugation is maintained without ^qfi- There are two special values of Gd which should be noticed. The 
first one is 1/46*0 = —^i,a{2m,0)\2=ii=A=o, with Gd > Go the stable diquark bound state can be found, and with 
Gd < Gq there is only unstable diquark resonance. The diquark self-energy in the normal phase can be separated 
into the vacuum part and the matter part as 

n(po,p) = n'"'^*(po,p) + nT=^=o(po + 2M,p) 

d^k l-(e + e0/2-/(gg')-/(e;+,) 



_ , [ ^ 

" 7 (27r)3 Po-?r-C^ 



r d^k 1 - e(-e') - e(-e) , , 

where the first term is a convergent matter term since the integrand function contains a Fermi distribution function 
for positive energy. To remove the divergence in the vacuum term, the self-energy is subtracted by n(2m, 0). Then 
one can defined a renormalized self-energy as W°'^{po,p) = li{po, p) — n(2m, 0). Meanwhile the coupling constant is 
renormalized as 

111 

(50) 



4Gr AGd 4Go' 

where Gr is related to the scattering length as \/AGr = m/47ra,, by taking the low energy limit p = k ^ m. The 
other limit for the coupling is Gc = — 11^ ^(0, 0)|t=/j=a=o- As Gd approaching Gc the mass of diquark bound state 
can become zero. If Gd > Gc the vacuum becomes unstable. Considering a system with total baryonic number 
density fixed, one obtain with the fermion number conservation, 

Ntot{t^,T) = NMF{f^,T) + NMSR{fi,T) 

= NMF{n,T) + \j ^ pJ^Y.^-^i^^^)M^) + Qi-^% (51) 

where Nmf is the fermion number density from free quarks, the total quark number density is given by Ntot = 
NcNf-^ with kp being the effective Fermi momentum. Together with the Thouless condition 

-^-ni,a(0,O)=0, (52) 
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Figure 8: Tc/Ep and iic/rn (left panel) and number density fraction (right panel) as as functions of normalized diquark coupling 
constant Gc/Gr- Here Tc and are the CSC critical temperature and chemical potential respectively. 



by tuning Go in the range [0, Go] for weak couplings and in the range [Go, Gc] for strong couplings, one can look at 
the BEC-BCS crossover along the CSC boundary {fic,Tc) with fixed kp. Because in the CSC boundary the gap is 
always zero, the six components of the diquark self-energy and the phase shift in above equations with indices i — R/I 
and a = 1,2,3 are the same. Then the diquark number density can be evaluated as 



NNSR{^J',T) 



p^dp 
'2^ 



duj ^^^^ P) [ns (^) + ^'(-'^)]- 



(53) 



When in strong couplings, Gr > 0, the diquark spectral density p(w, p) = Pci^, p) + Psii-^, p) consists of two parts, the 
continuous part (unbound part) and the pole part (bound part), with ps — Zb5{lj + 2^1 — Esp) + ZgS{uj + 2/i + Eg^) 
where ^^/b — ^McTITa^ I"=±Bb/bp-2m- ^B/Bp is the energy of the diquark and anti-diquark bound state with three 
momentum p. As a consequence the diquark number density is also divided into two parts, the unbound diquark 
contribution 7V„„ and the bound diquark Nb contribution. In the weak coupling case Gr < 0, only the unbound part 
is left. In Fig. [51 the parameters are set to A = 0.7 GeV, A^ = 0.65 GeV, m = 0.2A, and the total fermion number 
density is fixed at kp = 0.2m. With these parameters, one studies the diquark fluctuation effect near the unitary 
limit. The left panel shows the CSC critical temperature and chemical potential as functions of diquark normalized 
coupling constant. As the coupling becomes stronger, Tc rises and pc decreases. The effect of the diquark fluctuation 
is found to lower Tc and pc (solid curves) comparing to the mean field results Tmf and pmf (dotted curves). At 
a very weak coupling, the CSC boundary can be well approximated by the mean field approximation. The right 
panel shows fermion number density fraction of quarks and diquarks. The diquark component increases with the 
diquark coupling constant. In the strong coupling side with Gr > 0, Tc{pc) continuously increases (decreases) with 
the coupling constant. Both the free fermion fraction Nmf /Ntot and unbound diquark fraction Nun/Ntot decrease 
with the coupling constant once the effective chemical potential of the system becomes negative p — m < 0, while the 
diquark bound state contribution becomes nonzero and finally dominant. 

In this section, we summarize recent results on the diquark spectral densities in different regions of the phase 
diagram. The quark mass and CSC gap serve as boundaries for stable diquarks. In the 2SC phase, there are gap 
structures in the spectral density of the thermal diquark with red and green colors. The infinite peaks at (w, p) — (0, 0) 
GeV indicate five Goldstone modes. In the NJL model, the scattering length as of fermions can be determined by 
taking the non-relativistic limit for the T matrix, which is an advantage in comparison with the boson-fermion 
model beyond leading order calculation. Considering a total baryonic number density conserved system, tuning the 
scattering length along the CSC boundary where the Thouless condition is satisfied, the boundary of the CSC and 
the number density fraction of free quarks and diquarks are calculated. When the l/ag is negative and small, the 
diquark contribution to the number density is negligible and the system is in the pure BCS regime. The diquark 
unstable resonance has a remarkable contribution when the absolute value of the effective chemical potential is small. 
When l/us is positive, the diquark bound state will form. Once the effective chemical potential becomes negative, 
increasing 1 /as the number density fraction of diquark bound states becomes nonzero and quickly dominates and the 
system settles in a deep BEG regime. 
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V. BARYON AS QUARK-DIQUARK COLLECTIVE MODES 



Baryoii can be studied as a bound state of three constituent quarks |81l - l8a |. Quarks and gluons carrying color 
charges are confined inside baryons and mesons. The static properties of baryons in vacuum can be studied with the 
relativistic Faddeev equation [stI [88| . In this approach, the baryon is assumed to be stable with a separable form 
of the T- matrix. The original Faddeev equation can be reduced to the Bethe-Salpeter equation (BSE), which is an 
eigen-equation for the baryonic vertex. But with a nonzero temperature and baryon density, the Faddeev equation 
or the BSE is hard to solve numerically. The way around is to take the static approximation for the intermediate 
quark propagator in the Faddeev equation. Then the baryon can be studied at nonzero temperature by a two step 
process: first, the thermal diquark propagator is simulated by the Dyson-Schwinger equation, and second, the diquark 
is coupled with another quark to form a baryon. In previous works, the diquark is always assumed to be stable, 
but this is not necessarily true. Baryon can be a bound state of a quark and an unstable diquark, which is like a 
borromean state in nuclear and atom physics. 

To add the baryon field into the Lagrangian psp , a coupling term of the quark-quark-diquark-diquark is introduced 

as 

~ -^BB + ^BfB^^aVa^ + ^^a.^aPlj.B . (54) 

Here, tpafa — {''Pa'Pa) +f3a, and the baryonic field is defined as i? = Gb {ipa^a)- The terms of order 0(/3^) is neglected. 
Actually this is equivalent to take the static approximation in Faddeev equation. The baryonic fields in the NG basis 
are then denoted by B = {B,Bc)^ and B = {B,Bc). The baryon-quark-diquark vertices are Tbr = '^^ng and 

Tbi — *"72^3^''' respectively. The sum of the Lagrangians (I55t and ([51)) is the starting point for the further treatment. 

The 11-component in the NG space of the inverse baryon propagator is S^^ = —1/{2Gb) — S, where 

^(^) = --aT.I - K)[DR,a{K) + Di^,{K)] (55) 

is the 11-component of the baryon self-energy. The quark propagator in the NG space, Sfi, is diagonal in color 
space. In presence of a non- vanishing diquark condensate, S\i = S^i ^ Sfi- If the diquark condensate vanishes, 
Sli — — and Dji ^ = i'/.fc for any a, b. Inserting the spectral density form of the diquark full propagator Eq. 
(j44p into Eq. H55p . the summation over Matsubara frequency can be handled. Then the positive energy component of 
the baryon full propagator can be extracted with energy projectors, Sg^^{pQ,p 0) = ^Tr ['S'g^Ap^Q7*'] , where Ap 

is the energy projector Ap — ^[1 + s (707 • p + JqIMb) /Ep], with Ep — y/jP + M^ and s = ±1. In the homogeneous 
limit, p = 0, the energy projector assumes a simple form, Ap^Q — ^{1 + sjq), which is independent of Mb- The full 
expression of the positive energy component of the baryon propagator is. 



SB\ilo,0) = iTr{[-l/(2GB)-S>V} 



2 

= -1/(2G, 



fp'dp p duj 1 



ei4 + e;/(ei£^) + n(^) 



+ K,(.,p)-fp?(.,p)]«l±^}(l + e9i). (56) 

where summation is implied over e, ei. 

Fig. [9] shows the real and imaginary parts of the inverse retarded Greens function for baryons (positive energy 
component), again at points A,B,C, and D in the phase diagram of Fig. [S) The diquark coupling constant is taken 
to be weak, Gd = 3.11 GeV~^. The constant Gb = 10.04 GeV~^ is chosen to make the baryon mass of 0.94 GeV in 
vacuum. In the chiral symmetry broken phase with 7^ and A = (point A), there are no diquark condensates or 
resonances but stable baryon resonances: in the upper- left panel, one can see that ReS'^^(wB,0) = has a solution 
at lub + w 0.94 GeV, i.e., close to the nucleon rest mass. There is a region of ujb € [— 3(TOg -I- fiq),3{mq — /ig)] or 
Mb G [— 3mg, 3mg], where the imaginary part Im5^^(cjB, 0) is very small (smaller than 10~^ GeV) in the homogeneous 
limit. The position is just inside this region, i.e., Mb < Sniq: the baryon weighs less than its constituents. It is 
therefore stable, although its constituents by themselves are unbound, like in a Borromean state in atomic or nuclear 
physics. The upper-right panel shows the case with diquark resonances but outside the CSC phase (point B). There 



15 



W5 



0.8 
0.6 
0.4 
0.2 
0.0 
-0.2 
-0.4 
-0.6 
-0.8 



T=0.03,;i,=0.25 ( 




V 


\ 
1 
I 


-3-2-10 1 2 3 
iOB+3fig [GeV] 


T=0.15,;i,=0.36 ( 




\ ^ 

\ 





-2-10 1 2 
wb+3/1, [GeV] 



0.8 
0.6 
0.4 
0.2 
0.0 
-0.2 
-0.4 
-0.6 
-0.8 



1.0 
0.8 
0.6 
0.4 
0.2 
0.0 
-0.2 
-0.4 
-0.6 
-0.8 
-1.0 



T=0.03,;i,=0.33 ( 


JeV 


\ 
\ 

\ 


\ 
\ 


-3-2-10 1 2 3 
wb+3//, [GeV] 


T=0.03,/x,=0.36 ( 

^ \ ^ 


JeV 




\ 



-2-10 1 
wb+3//, [GeV] 



Figure 9: The real (blue dashed) and imaginary (red solid) parts of the inverse propagators for baryons as functions of energy 
Lo at different T and fiq. From top to bottom, the first panel: T = 0.03 and Hq — 0.25 (point A). The second panel: T — 0.03 
and Hq — 0.33 (point B). The third panel: T = 0.15 and fiq = 0.36 (point D). The fourth panel: T — 0.03 and fiq = 0.36 (point 
C). All units in GeV. 



is no positive energy baryon pole in this case. In the region of higher temperatures and quark chemical potentials 
where chiral symmetry is restored and where there are neither diquark condensates nor resonances (point D), there 
are also no baryon resonances and the absolute value of ImS'^^^ is very large. This case is shown in the third panel. 
In the CSC phase (point C), there are baryon poles but with nonzero imaginary parts, indicating unstable baryon 
resonances, as shown in the fourth panel. 

There are a lot of works using the simplified Faddeev equation by static approximation to study baryon properties in 
vacuum and nuclear matter, see, for example, Ref. (90l[9]|. But the issues the authors focus on are different. To give 
the diquark propagator they used the proper time regularization method which introduces an effective confinement, 
but the method is not applicable to nonzero temperature case. The diquark T-matrix is approximated by the constant 
term plus the pole terms, which is equivalent to the assumption of a stable form diquark propagator, while as 
presented in Fig. [5J the baryon can also be formed by an unstable diquark and a quark. The pole approximation for 
the diquark propagator will miss some of important physics like Borromean state. In further calculation the physics 
mass of baryon bound state or baryon resonance can be obtained at any given T and fi on the phase diagram based on 
the mean field approximation by the NJL. In the present approach, at low temperatures, the baryon mass have only a 
slight decrease when the chemical potential rises, but in Ref. (ool . [9l| baryon mass decreases significantly. The reason 
is that the vector meson is not included and a large baryon number density could not be obtained by increasing the 
chemical potential. Also the confinement mechanism at finite temperatures should be incorporated. In Ref. [92\ , the 
static approximation and the stable diquark are used. The authors considered a three-flavor NJL model, where the 
baryon mass is found to decrease by 25% at normal nuclear matter densities. These issues can also be considered in 
the present framework. 

In summary, diquark propagating modes are derived with an NJL-type model in different regions of the phase 
diagram of strongly interacting matter. Baryon formation and dissociation in dense nuclear and quark matter is then 
studied via the baryon poles and spectral densities, incorporating the previously obtained diquark propagator. The 
stable baryon resonances with zero width are present in the phase of broken chiral symmetry, where the diquarks 
could be an unstable resonance. This indicates that the baryon can be a Borromean like bound state. There are no 
baryon poles in the chirally symmetric phase. In the CSC phase, baryon poles exist, but they are found to be unstable 
due to a sizable width. 
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VI. NONLOCAL EXTENSION OF NJL MODEL AND ITS APPLICATIONS 

The NJL-type model applied to quarks is a successful schematic effective theory for QCD, in which the interaction 
between quarks are described by point-like couplings. The model can be used to study the spontaneous chiral symmetry 
breaking and quark pairings. But the basic version of the NJL model have some shortcomings. For example, the 
NJL model can not exhibit confinement as in QCD. The other one is the constituent quark mass is independent of 
momentum which is in conflict with the lattice data and the Dyson-Schwinger equation from QCD. The second point 
can be studied by a nonlocal extension of the point-like coupling in the NJL model. The nonlocal quark current reads, 

Ja/ = jd'^xid!^X2f{xi)f{x2)q{x-xi)TMq{x + X2), (57) 

in which = cr, vr stand for the scalar meson and the pion. T„/^ are 1 and 7^ra respectively with flavor indices 
a = 1, 2, 3. / is a form factor in momentum space, usually taken as a Gaussian type /^(p^) — exp(— p^/A^). The basic 
NJL model is non-renormalizable and some form of ultraviolet regularization is necessary with a cutoff parameter 
which is a part of the model. But in the nonlocal NJL model the momentum integration are automatically convergent 
in the loop diagrams and no additional regularization method is needed (93i] . 

The parameters of the model including the coupling constant G, the quark current mass mc and the momentum 
cutoff A are determined by fitting the pion mass Mt^ and the decay constant in vacuum. The meson sectors are 
defined by introducing the auxiliary scalar field a — G {Ja{x)) and pseudo-scalar Tr°' — G (J^(x)). The dressed quark 
propagator in the mean field approximation is determined by the following equation 

S{p)-' = f-m,~^{p), (58) 

where S(p) is the self-energy which turn out to be S(p) = nidf'^ip'^), with rrid = iGTM Tr[rj\/5(fc)/^(fc^)] a 
momentum independent constant serving as an order parameter for the dynamical chiral phase transition. The 
IPI diagram of the self-energy vanishes due to the integration for xi,X2 goes from —oo to +oo. Comparing to 
the classic NJL model, the constituent quark mass now depends on the three momentum by a Gaussian factor. In 
Dyson-Schwinger equation from QCD, the dressed quark propagator has the form S{p) — Z{p'^)[f — Af(p^)]. The 
renormalization function Z{p'^) can also be obtained in the nonlocal NJL framework by considering the thermal meson 
correlation beyond mean field approximation [93]. 

The meson propagator is given by RPA for the Bethe-Salpeter equation D^^^^[p) = —G^^ + Il^\p), with 11*^ (p) 
polarization function in the mean field approximation |93|. 

U"'ip)=t [ f[{k+p/2)^]f[{k-p/2f]Ti[S{k+p/2)rMS{k-p/2)rM]. (59) 
Jk 

The pion mass is obtained by the pole condition — G^^ + n''(p^ = mj) = 0. Then the quasi-particle propagator can 

be written as D^'ip) = ^33^ with = gn'(p)/Op^ \p'^=ml- 

To calculate the pion weak decay constant, the weak current is introduced by a delocalization procedure for the 
quark fields [9^ , that is 

qiy)^Q{x,y)^Eix,y)q{y), (60) 

where E{x,y) — Vexpi dz^^[V^{z) +y^°(z)7^]T° is the Schwinger phase factor, V^{z) and A'^{z) are vector and 
axial-vector gauge fields respectively. The nonlocal current is modified as 

Jm = Jd*xid'^X2f{xi)f{x2)Q{x-xi,x)TMQ{x,x + X2)- (61) 

Then the weak vertices are introduced in the present nonlocal NJL model of two types: the weak current coupled with 
a quark and the weak current coupled with the quark meson vertex r^^^. Both give rise to the bubble diagrams 
which contribute to pion weak decay. 

The nonlocal extension of the quark NJL model is inspired and therefore refiects some important features of QCD. 
The momentum dependence of the constituent quark mass can be considered and the original sharp momentum cutoff 
is replaced by a smooth one. Then all the loop diagrams are automatically convergent and hence the cutoff dependence 
of the model is highly weaken. The last point can be looked at from the calculation of the meson loop effect in the 
1/iVc expansion. In the nonlocal study [93, the next to leading order contribution to the quark condensate turn out 
to be positive in vacuum, in contrast to the result obtained from the local NJL model [94] , in which the beyond mean 
field correction is a nonlinear function of the momentum cutoff Am- 
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VII. GINZBURG-LANDAU EFFECTIVE THEORY IN THREE-FLAVOR DENSE QUARK MATTER 

WITH AXIAL ANOMALY 

The QCD phases in high density regime are controlled by the chiral and diquark condensates = (qq) and d = (qq). 
The interplay between the Nambu-Goldstone (NG) and color superconducting (CSC) phases can be described in an 
model-independent way by the Ginzburg-Landau (GL) theory (osl . [o^ . It has been predicted in the GL theory that 
a new critical point and smooth crossover arise from the coupling between chiral and diquark condensates induced 
by axial anomaly at the low temperature in the QCD phas e diagram [97l - [99| . Such a coupling is also related to 



the continuity between the quark and hadroni c matter |10C| . The new critical point can also be confirmed in the 
three-flavor N JL model with axial anomaly [tqI . llOll Il02l | . 

The form of the GL free energy to the sixth order in the fields can be constrained by the symmetry, 

g ^ SU{3)l <E) SU{3)b. <E) U{1)b <E) U{1)a <E) SU{3)c (62) 

where subscripts L, R, B, A and C stand for left-handed flavor, right-handed flavor, baryon number, axial charge, 
and color symmetry. The left-handed and right-handed quark fields transformed under Q as 

qx ^ e-'^'^e-^'^VxVcqx, (63) 

where A = L{+), R{—), Vi^/jf/c ^i'^ rotational matrices of SU{3)[^/f;/c, and Oj^/g are rotational angles of U{1)a/b- 
The chiral fields are defined and transformed as 

<j> ^ e-^'^^VL^Vl (64) 

where i,j denote the flavor indices and a denotes the color indices. We see that (f> is invariant under Z{2)a C U{1)a- 
Then we have transformation rules for these field quantities, 

_^ VLcf><l>^Vl 

detcj) e-^^^-*det0. (65) 

We see that det(j> is invariant under Z{6)a C U{1)a- 
The diquark fields are defined as 

^ '^"/97e^J'c ((gcA)^i (^a)^) = ea/37e»ifc (('?^)^C((7A)^;) , (66) 

where C — «7^7o- Under Q, the diquark fields transform as 

dx e'^'^'^^e^^'^^VxdxV^, 

d{ ^ e'^'^'^^e~^'^^VSd{vl (67) 
The color singlet quantities in diquark fields transform as 

dxdl ^ Vxdxd{vl 
dxdl^ ^ e^'^o^Vxdxdl^vl^, 

detdx e'^'^'^^e^'^'' detdx- (68) 
Then the most general form of the GL free energy which is invariant under the transformation of Q read 
Vt{(j),dL,dR) = VL^((I)) + Vld{dL,dR) + Vl^d{<i>,dL,dB), 

n^{cf>) = ^Tr</.t0+|(Tr</>t0)2 + |Tr(0t0)2_|(det0-Kdet</.t) 

nd{dL,dR) = aoTr[di4 + d^^] + /?i([Tr(did[)]2 + [Tr(dfl4)]2) 
+l32[TiidLdlf + Tr(dfl,4)2] + /33Tr[(dfl4)(di4)] 
-f/34Tr(di4)Tr(dfl4), 
n^d{c^,dL,dB,) = jiTr{dBd[c^ + dLdy^) + XiTr{dLd[c^(f>^ +dndy^(j)) 

+X2Ti{dLdl + di?4) ■ Tr(</'t0) + A3[det<?!> • Tr(dL4'/'^) + ^.c.]. (69) 
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For three massless flavor case, the most symmetric condensate is in the form 



Then the GL free energy reads, 



4) = diag(cr,cr,CT), 
dL = -dfi = diag(d,d,d). (70) 



= -^a ^-a +-a + -a 

+ ^d' + ^d*-j<fa + \cfa^ (71) 



where coefficients a,b,c, f,a(3^X come from those in Eq. (|69p. Here a, a are two essential parameters to drive the 
phase transition, b can change sign with T, /x, so a positive term (/ > 0) is introduced to stabilize the system 
for negative b. /? is positive definite from effective theories and weak-couphng QCD. The and d^cr terms are from 
axial anomaly so the coefficients c and 7 are related and are all positive. A is also positive from the NJL model and 
weak-coupling QCD. 

There are four phases: the normal (NOR) phase with a ~ and d ~ 0, the CSC phase with ct = and d ^ 0, the 
Nambu-Goldstone (NG) phase with cr ^ and d = 0, and the coexistence (COE) phase with cr 7^ and d ^ 0. The 
phases at a specific set of parameters can be determined by comparing the global minima of the free energies flpfORj 
^csc, ^NG and ftcoE- From NOR or CSC phase with ct = to the NG phase with ct ^ 0, there is a first-order 
transition phase for the chiral symmetry breaking or restoration. The transition between the NOR phase with d = 
and the CSC phase d 7^ is a second-order one with a discontinuity of 

VIII. NJL MODEL FOR THREE-FLAVOR DENSE QUARK MATTER WITH AXIAL ANOMALY 

The effective Lagrangian for three-flavor dense quark matter with axial anomaly can also be derived from the NJL 
model (tqI . [soj . The NJL Lagrangian reads 

£ = q(^7^9^-m,+//7°)g + /:(4)+£(6)^ (72) 
where — + >C^^''. The four-fermion interaction term reads 

^cf' ^ ^ X! [iqil5TA>^A'qc){qcil5TA>^A'q) + iqTA>^A'qc){qcTA>^A'q)] , (73) 

A,A'=2,5,7 

where qc — Cq^ and qc = q^C with C = «7^7o- The six-fermion interaction term is C^^^ = + ^cl^j^jj^g, where 
is the standard 't Hooft term, 

4'^ = -i^(det[g(l-H75)9]+ det[g(l-75)g]| 

k-ilaovr navor ) 

= -Ksijk {[gi(i + 75)%]fe(i + 75)gj][g3(i + 75)%] 

+ [91(1 - l^)q^][q2{l - 75)gj][g3(l - 75)'Zfe]} , (74) 

where 91,2.3 denote u, d and s respectively. Here we introduce a mixing term C^^^^^^^ for the coupling between the 
chiral and diquark condensates, 

4ting = \k' J2 {i^PBA{l + l,)q][qrB\A'il+l,)q''][fr^\A'{l + l,)q] 
A,B,A'=2,bJ 

+ [qTAXA'{l - l5)q'^][fTB\A'{l - lz)q][qPAB{l - 75)9]} , (75) 



where Pab are matrices in flavor space, 
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Si life we can decompose the chiral and diquark condensates in the following form 

^ = -qq= A,) ® (.^f ® (0f " + 0f "75 + "7'^ + 0;f5"7''75 + "^^^^ 5''^) 

D = -qf = (C'°'-A„) ® (4-°V/) ® (rff " + (f5P'"75 + d^r^^ + d'^ti'^^r, + tf^P-^"^) 

Then we have 

qTfq = -Tr[gfgr/] = 2(?!)^, 

975T/g = -Tr[gg75T/] = 2</.^, 

q'^i'YsTAXA'q = -Tr [gg^nsr^AA'] = 2irff ^, 

9«75TaAa'«'^ = («'^«75TaAa'«)^ = -2i(rff 

g^TAA^'? = -TV [qfrAXA'] = 2df^, 

qTAXA'q'' = {(pTA\A'q)^ = 2{d^'^)^. 

In the mean-field approximation, considering the CFL channel and dropping the 0~ state, we have 

{cj,i)Tf = a, (05^)r/ = 0, 

= \dbA'A, {d^'^) = 0. 

Now the Lagrangian becomes 

^ 4Gaqq-6G(7^, 
Cf ^ H [d*{f^5TAXAq) - {qi5TA\Aq'^)d] - 3^f|rf|^ 
£(6) _^ -2K(j'^qq + AKa^, 

4ting ^ -^Ml'a? - [d*{fTA\Al,q) - {qTA\Al,q'')d] + 



Thus the Lagrangian can be expressed in the Nambu-Gorkov basis \1/ = ^^(9, ?^)^, 



2 



where 



with 



[/ = 6G(j2 + 3i?|d|2-4i^a3-^X'a|dp. 



M = mq - 4G + 2/ra2 + 
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Table I: The parameters for the gap equations without 't Hooft term used in Ref. (79l |. 





mq [MeV] 






Kh:' 


M [MeV] 


cr^/^ [MeV] 


I 





1.926 


1.74 


12.36 


355.2 


-240.4 


II 


5.5 


1.918 


1.74 


12.36 


367.6 


-241.9 



Thus we obtain the thermodynamic potential as 



VL = - 



■E 



Wo 



16Tln l-fexp(-^) 



2Tln(l + exp(-^) , 



T 

+ U{a, d). 



(84) 



where = ^/{Ep ± + |Ap, ujf = ^(£'p ± -f 4|A|2 with = + K'P{(t, d). Thus the gap equations are 



^ + [-2/K=') + 1] 



dd 



+ [-2/K±) + 1] 



da 

dujf 
~dd 



da ' 



9d 



= 0. 



(85) 



The gap equations without mixing 't Hooft term has been solved with parameters in Table |T] [79] . There are three 
phases in the phase diagram: the NOR phase with cr = 0, d = 0; the NG phase with cr 7^ 0, d = 0; the CSC phase 
with cr = 0, d 7^ 0. If quarks are massive 7^ 0, there is a critical point at high temperature which is called the 
Asakawa-Yazaki point, which is the endpoint of the first-order phase transition line in the phase diagram leading to 
a crossover. 

When the mixing term is introduced with a strong enough chiral-diquark interplay K' = 4:.2K, a new critical point 
at low temperature will emerge just as predicted in the Ginzburg-Landau approach [97l[98|. The interaction between 
chiral and diquark condensates weakens the chiral symmetry spontaneous breaking and leads to the COE phase with 
CT 7^ 0, d 7^ at low temperature. As a consequence, the first-order phase transition between the NG and CSC phase 
becomes a crossover. This is a new critical endpoint at the other end of the first-order phase transition line. 

At the same time, the mixing term also induces a BEC-BCS crossover (ssl - is^ [6l|. Here a new criterion in the 
dispersion relation is used to define a BEG state. For /i > M, there will be non-vanishing momentum p = \/ [i^ — M'^ 
to give the minimum energy and it is a BCS state. But for /i < M, the minimum energy has to at p = which means 
the system is in a BEG state. Since we have the COE phase with non-vanishing chiral condensate, the BEG state 
occurs in phase diagram on the left side of the curve /x = M{^, T). Physically, the BEG state can also be explained as 
a compound particle including two strong coupling quarks because the diquark channel is strengthened by the mixing 
term just as H' = H + \K'a for sufficiently large K' . 

However, it was pointed out that the 2SC is present if K' is sufficiently strong [sO*]. This is due to that the axial 
anomaly induces a mutual amplification of the str ang e chiral condensate and the non-strange diquark condensate. As 
the consequence, the critical point found in Ref. [79| only survives for very narrow parameter space of K' otherwise 
most parts are covered by the 2SC phase fso!]. 



IX. GINZBURG-LANDAU APPROACH TO COLLECTIVE MODES IN SPIN-ONE CSC 

The Ginzburg-Landau approach can be used to study the collective modes in spin-one color superconductors (sgI . 
IstI. |43|. The spin-one color superconductor involves pairing of quarks of same flavor. The diquark condensate or the 
order parameter A is then a color anti-triplet and spin triplet, so it is a 3 x 3 complex matrix and transform as 



A ^ J7Ai? (86) 

where U — exp(i6'(iA(i) G U(3) = SU(3)cXU(1)b and R — exp{iaiJi) e 50(3) are transformation matrices. Here are 
eight Gell-Mann matrices and Aq is normalized unit matrix. {Ji)jk = —i^ijk are generators of S0(3)pi, 9a{a = 0, 8) 
and ai{i = 1,2,3) are rotation angles in U(3)l and S0(3)r group space. There are 18 real parameters in A, among 
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which 12 parameters are carried by the U(3)l x S0(3)r transformation making a 12-dimensional degenerate vacuum 
manifold. Then Acan be parametrized by the remaining 6 real parameters which characterize different vacuum states, 

Ai iSs —iS2 

A = I -iSs A2 i<5i I . (87) 
(§2 -iSi A3 

A. Ginzburg-Landau free energy and ground states 

Up to fourth order in A and two derivatives, the most general U(3)l x S0(3)r and parity invariant Ginzburg-Landau 
free energy density functional can be written as 

J^[A] = aiTr(a,Aa,At) + a2{d,Aai){djAl^) + 6Tr(AAt) 

+di[Tr(AAt)]2 + d2TV(AAtAAt) + rf3Tr[AAT(AAT)t]. (88) 

The time-dependent GL functional, or Lagrangian, is then written as 

£ = iciTr[AtaoA] -f c2TrPoAt)(aoA)] - J-[A], (89) 

The ground state is found by minimizing -^[A]. The phase structure, or orientation in the field space, of Ag depends 
on d2 and ds. Here, only the following four ground states represented by the following four matrices occupy a part of 
the phase diagram. 



^poloar 





A, = /3 , (90) 



where a = {d2 + d3)/[2{2d2 + da)] and /3 = y/d2/{2d2 + da). The pattern of spontaneous symmetry breaking deter- 
mines the low-energy spectrum of the system, i.e., the NG bosons. While some of the NG bosons are associated with 
the generators of the color SU(3)c group and arc thus eventually absorbed in gluons via the Higgs- Anderson mech- 
anism, those stemming from spontaneous breaking of baryon number or rotation symmetry remain in the spectrum 
as physical soft modes. As we will now see, some of the phases exhibit the unusual type-II NG bosons, in accordance 
with general properties of spontaneously broken symmetries in quantum many-body systems. 

B. CSL phase 

When d2 + ds > and d2 > d^, the ground state is the CSL phase in which the spin and color are coupled in the 
pairing, so the symmetry breaking pattern is U(3)l x S0(3)r — ^ S0(3)v. There are 9 broken generators leading to 9 
NG bosons as follows, 

• Ao «) 1. Type-I NG singlet, E"^ ~ (ai + a2)P. 

• V^('^7 (g) 1 - 1 ® Ji), ^^(-As (g 1 - 1 (g) J2), ig) 1 - 1 (g Js)- Type-I NG triplet, E"^ ~ (ai + 02)^2. 

• Aa 1, a = 1, 3, 4, 6, 8. Type-I NG 5-plet, E'^ ~ (ai + 02)^^. 

C. Polar phase 



When ds < and d2 + d^ < 0, the ground state is the polar phase. The symmetry breaking pattern is U(3)l x 
S0(3)r U(2)l X S0(2)r. There are 7 broken generators which, however, give rise only to 5 NG bosons, organized 
in the following multiplets. 
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• ® 1, where 7^3 = "^(-^o — V^Ag) = diag(0,0, 1) is the projector onto the third color. Type-I NG singlet, 

£'2 ~ aifci + (ai + a2)A;|. 

• 1 (g) Jj, j = 1, 2. Type-I NG doublet, ~ (ai + a2)k\ + oiA;^. 

• Ao, 1, a = 4, 5, 6, 7. Type-II NG doublet, E'^ ~ o^fci + (ai + a2)^fc|. 

The presence of type-II NG bosons is due to nonzero color density of the polar ground state. 

D. A-phase 

When (^3 > and ^2 < 0, the ground state is the A-phase. The symmetry breaking pattern is U(3)l x S0(3)r — )• 
U(2)l X S0(2)v- Unlike in the polar phase, the diquark spin is now circularly polarized. Among 7 broken generators, 
there is only one giving rise to a type-I NG mode, 

• y^(^3 ® 1 + 1 ® Js)- Type-I NG singlet, E"^ ~ (ai -h a2)h\ -h aifc|. 

The rest 6 generators produce only 3 type-II NG bosons due to non-zero color and spin density of the A-phase vacuum, 

• Ao; O 1, a = 4, 5, 6, 7. Type-II NG double, E"^ ~ (ai -h aifk\^ a\k\. 

• 1 \{Jx ± U2). Type-II NG singlet, E"^ ~ af yfci (ai a-ifk\. 



E. e-phcise 

When da > (i2 > 0, the ground state is the e-phase. The symmetry breaking pattern is U(3)l x S0(3)h — > 
U(1)l X S0(2)v. The spin of the second diquark color is longitudinal polarized, while that of third color is circularly 
polarized. Out of the 10 broken generators only two correspond to type-I NG modes: 

• \plVi ® 1, where Vi = diag(0, 1,0) is the projector onto the second color. Type-I NG singlet, E'^ ^ aifcj^ -|- 
(ai + a2)A;|. 

• \/i(^3 O 1 -I- 1 (g) J3). Type-I NG singlet, E'^ ~ (ai + 02)^5. + aikj. 

The remaining 8 generators give rise to 4 type-II NG modes due to nonzero color and spin density of the e vacuum, 

• ^(Ai ± iX2) (g) 1. Type-II NG singlet, E'^ ~ afk^ + (ai + a2)2fc|. 

• ^(A4,6 ± «A5,7) (g 1. Type-II NG singlet, E"^ ~ (ai + a2fkl + afk^ 

• 1 (g) -l=(Ji ± i J2). Type-II NG singlet, E"^ ~ ajkl + (ai + a2fkl 

In summary the low-energy physics of spin-one CSC is analyzed in terms of the NG excitations within the Ginzburg- 
Landau theory. The four phases that appear in the phase diagram possess different NG modes of the spontaneously 
broken color, baryon number and rotational symmetry. Those stemming from the color symmetry will eventually be 
absorbed into gluons, making them massive by the Anderson-Higgs mechanism. The other NG bosons will remain 
in the spectrum as physical soft modes. Unlike in all the other phases, in the isotropic CSL phase all quarks can be 
gapped so that the NG bosons are the only truly gapless states in the spectrum. 

X. SUMMARY 

We give an overview on recent progress in quark pairings in dense quark matter. These progress include the 

BCS-BEC crossover in the CSC with and without external magnetic field, baryon formation and dissociation in 
quark/nuclear matter, Ginzburg-Landau effective theory on dense quark/nuclear matter with anomaly, and collective 
and Nambu-Goldstone modes in spin-one CSC. 

The boson-fermion model in the cold atom system can be extended to relativistic case to describe the relativistic 
BEC-BCS crossover in the CSC. In a charge conserved system, by tuning the bare boson mass the effective coupling 
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between fermions changes, and the fermion number fraction and fermionic chemical potential as functions of the 
effective coupling also change indicating a crossover between BCS and BEC states. The pseudo-gap and effects of 
thermal bosons can also be systematically studied within the CJT formalism. The BCS-BEC crossover can also be 
dealt with in a pure fermionic model in which the fermion scattering length can be derived in comparison with the 
result at the low momentum limit for the T-matrix. The bosonic degrees of freedom are introduced by a bosonization 
procedure. In dense quark matter with fixed total baryonic number density, tuning the coupling constant along the 
CSC boundary, the baryon number fraction of free quarks and thermal diquarks vary with the scattering length, 
indicating a typical BCS-BEC crossover. Both the double channel (boson-fermion model) and single channel (pure 
fermionic model) models can describe BCS-BEC crossover in relativistic case equally well. The strong magnetic fields 
which exist in pulsars and non-central heavy ion collisions can tune the BCS-BEC crossover in a charged fermion 
system, but the origin is different from that in a cold atom system where the magnetic field is used to change the 
fermion coupling via Feshbach resonance. A strong magnetic field up to the energy scale of QCD have impact on the 
fermion energy dispersion. Sitting in a BEC state, varying the magnetic field can induce an oscillation of the fermion 
number fractions which makes the system be in the BCS- and BEC-dominant region. As the magnetic field becomes 
even stronger all fermions are pressed to the lowest Landau-level, hence the system settles down in the BCS regime. 

Baryons can also be regarded as a bound state of quark-diquark coupling in the NJL-type model, which is a 
simplification for the Faddeev equation. The diquark spectral density can be obtained by solving the DSE, with which 
the full diquark propagator can be written in a spectral density form and the summation over Matsubara frequencies 
can be analytically done. Then coupling the diquark with a quark within DSE, the baryon spectral density can be 
calculated in different phases. The formation and dissociation properties of baryon can then be investigated. In 
previous studies, diquarks are always assumed to be stable or as quasi-particles. This is not necessarily true because 
a baryon can also be a stable bound state of a quark and an unstable diquark, which bears some similarities to 
Borromean state in nuclear physics. The saturation of nuclear matter is still a challenge in current model. First one 
should include the vector meson channel to obtain a large enough baryonic number density. On the other hand, an 
effective confinement should also be included. In the normal quark NJL model, the absence of confinement is obvious 
since the quark mass is momentum independent. The nonlocal extension of the NJL model may provide a mechanism 
to include the confinement, in which the constituent quark mass is a momentum function. The dynamic quark mass 
and renormalization factor calculated within the nonlocal NJL model have a great consistency with the DSE analysis 
and the lattice results. The other advantage of the extended model is that in the bubble diagrams there always exist 
a Gaussian form factor which makes the integrals automatically converge. 

The phase diagram at low temperatures near the confinement and CSC boundary is still not clear. Recently with 
the three-flavor NJL model including the axial anomaly, a low temperature critical point near the baryon chemical 
potential axis is found in a Ginzburg-Landau analysis as a result of the interplay between the chiral and diquark 
condensates. One can derive the Ginzburg-Landau effective potential with the axial anomaly term within the NJL 
model by introducing the six-fermion mixing term of the diquark-diquark-quark-anti-quark coupling. When coupling 
constant of the mixing term is strong enough, a new critical point at low temperature can emerge just as predicted 
in the Ginzburg-Landau approach. But it has been argued that the 2SC is present instead of the critical point if the 
coupling of the mixing term is strong. This is due to that the axial anomaly induces a mutual amplification of the 
strange chiral condensate and the non-strange diquark condensate. 

The low-energy physics of spin-one CSC is analyzed in terms of the NG excitations within the Ginzburg-Landau 
effective theory. The NG modes of the spontaneously broken color, baryon number and rotational symmetry are 
analyzed for four typical phases. Those stemming from the color symmetry will eventually be absorbed into gluons, 
making them massive by the Anderson-Higgs mechanism. The other NG bosons will remain in the spectrum as 
physical soft modes. Unlike in all the other phases, in the isotropic CSL phase all quarks can be gapped so that the 
NG bosons are the only truly gapless states in the spectrum. 
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